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ABSTRACT 



A two-dimensional, two-passage simulation of the relative flow through a transonic 
fan at M = 1.4 was designed with a view to providing an apparatus in which to assess the 
effectiveness of passive vortex generator techniques in alleviating shock-boundary layer 
interaction effects. The design of the model and the results of six tests in the transonic 
cascade blowdown wind tunnel are described. Schlieren photographs of the shock 
structure were obtained at back pressures lower than the design value. The back-pressure 
control valve was identified as being critical to completing the experimental simulation. 
The flow through the cascade geometry was computed at design pressure ratio using an 
Euler code. Modifications to the grid are recommended before thin-layer Navier-Stokes 
calculations are performed. 
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I. INTRODUCTION 



A. SHOCK-BOUNDARY LAYER INTERACTION 

The motivation for the present study is as introduced in Ref. 1: 

Outboard of approximately 50-60% span on current fan designs, the peak relative 
Mach Number reaches 1.3. There is a normal shock wave that stands in front of the 
passage that is of sufficient strength to separate the boundary layer on the suction side, 
as shown in Fig. 1. Depending on the strength of the shock-boundary layer interaction, 
the boundary layer may or may not reattach itself to the suction side. As a result of this 
separation, there is little or no pressure recovery downstream of the shock and there are 
high total pressure losses. If this separation is able to be suppressed, which would 
improve the boundary layer characteristics downstream of the shock, the stage pressure 
rise would increase and the fan efficiency would improve. In addition, the fan rotation 
speed could increase, allowing a reduction in size and weight and more efficient turbine 
speeds. 

A diagram of the basic flow features of the shock-boundary layer interaction are 
shown at the bottom of Fig. 1. The step increase across the shock cannot be handled by 
the subsonic portion of the boundary layer, therefore, this pressure rise is partially 
transmitted upstream through the subsonic region. This results in the divergence of the 
stream lines and formation of compression waves upstream of the shock. If the shock 
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pressure rise is sufficiently strong, the boundary layer will separate and cause a more 
severe thickening of the upstream boundary layer. This will cause the coalescence of the 
compression waves into the leading shock. In order for the flow to turn back towards the 
axial direction, a rear shock is formed. The leading and rear shock combination is called 
the lambda foot. The total pressure loss through the lambda foot is lower than the loss 
across the normal shock, resulting in a free shear layer originating at the shock triple 
point. The undisturbed boundary layer thickness at the shock location, 8 0 , is the 
normalizing length used in describing the geometry of the interaction. 

Recent experiments have examined several concepts which have been proposed for 
suppressing the induced separation caused by the shock-boundary layer interaction. Three 
methods investigated at the United Technologies Research Center (UTRC) include low- 
profile vortex generators, a porous wall with a passive cavity and surface contouring [Ref. 
1]. Additional methods examined at the NASA Langley Research Center include large 
eddy breakup Devices (LEBU) at small angle of attack, spanwise cylinders and vortex 
generator jets [Ref. 2]. The most promising and adaptable for modem transonic cascades 
are the vortex generators and the vortex generator jets. 

B. VORTEX GENERATING DEVICES 

One proposed method for supressing or eliminating the shock-induced separation 
is by the replenishment of the boundary layer momentum near the surface [Ref. 2]. 
Replenishment can be effected by a redirection of the momentum in the free stream or 
outer region of the boundary layer into the near-wall flow. A mechanism for effecting 
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the momentum transfer is via the introduction of streamwise vorticity. Such a momentum 
exchange in the boundary layer can enable the layer to adjust to the sudden pressure rise 
without separation [Ref. 1]. Various techniques can be used to introduce streamwise 
vorticity. Two techniques, vortex generators and vortex generator jets, are described in 
the following sections. 

1. Vortex Generators 

Passive vortex generators with a height on the order of the boundary layer thickness 
have long been known to increase the mixing between external streams and the boundary 
layer. Such vortex generators, producing relatively large streamwise trailing vortices, 
have been used effectively to delay boundary separation, to enhance aircraft wing lift and 
to avoid or delay separation in subsonic diffusers [Ref. 2]. However, because of the 
fullness of the mean velocity profile of high Reynolds Number turbulent boundary layers, 
a reduction of the height of the vortex generators to a fraction (25-50%) of the boundary 
layer height has been proposed [Ref. 2], The small profile offers smaller perturbations 
to the supersonic freestream. Also, they are rugged in design, capable of withstanding 
centrifugal loads and other harsh conditions. Lastly, they have demonstrated a working 
flow incidence of ±5° [Ref. 1], The low-profile vortex generators, shown in Fig. 2, were 
tested in a flow in which the boundary layer was approximately 1.28 in. thick. They 
include the doublet and the wishbone [Ref. 2]. Other designs include the asymmetric 
candleflame and the singlet (a doublet with the second wedge removed) [Ref. 1]. 

UTRC investigated the effects of low-profile vortex generators on suppressing 
shock-induced separation in a slowly diverging circular duct. Details of the experiment 
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are found in Ref. 1. The shock position within the duct was controlled using a translating 
plug valve following the test section. Wall pressure distributions through the interaction 
region were obtained and surface flow structure was observed using an oil flow technique. 

The results showed that the pressure rise increased towards the ideal value when 
vortex generators were attached upstream, and the region of the pressure rise was 
shortened. Thus the shock-induced separation was substantially reduced. Depending on 
the configuration of vortex generators, the shock-boundary layer interaction length could 
be shortened by approximately 50%. Also, the displacement and momentum boundary 
layer thicknesses downstream were substantially reduced. A sketch which illustrates 
changes in the pressure recovery obtained with vortex generators is shown in Fig. 3 in 
comparison with the ideal pressure recovery and the baseline pressure recovery. 

Of the remaining methods of suppressing shock-induced separation investigated in 
Ref. 1, the porous wall with the passive cavity did reduce the shock loss, but the 
boundary layer development was aggravated. Also, significant mixing losses occured 
downstream. The use of surface contouring moderately reduced the boundary layer 
thickness, but did not affect the size of the separation region. 

2. Vortex Generator Jets 

The use of vortex generator jets was first proposed by Wallis [Ref. 3] in Australia. 
The vortex jets are usually, but not as intended here, an active method where jets of air 
are blown through holes in the solid surface. The holes are inclined at an angle to the 
surface and skewed with respect to the freestream. A schematic of the arrangement 
(taken from Ref. 2) is shown in Fig. 4. The arrangement of the jets on the surface is 
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much like the arrangement of solid vortex generators. The relatively weak counter- 
rotating vortices that are formed within the jet fluid in crossflow are replaced by a 
stronger streamwise vortex trail in the ambient fluid downstream of injection and close 
to the surface [Ref. 2). 

Vortex generator jets have proven effective in reducing the extent of the turbulent 
boundary layer separation resulting from adverse pressure gradients. However, the status 
is summarized in the following points stated in Ref. 4: 

• The feasibility of the technology is not firmly established 

• It is possible to adapt the arrangement passively to fans because of the cascade 
arrangment 

• There is a negligible drag penalty of the vortex generator jets when compared to the 
vortex generators when the jets are off. 

The effect of skewed jets and their relative angles (0) to the freestream on 
separation in a subsonic diffusing flow were investigated in Ref. 4. The jets were skewed 
45° to the vertical and 0 was varied from 0° (downstream) to 180° (upstream). The 
following conclusions were reached. 

• Substantial reduction of the size of the stall region size occurred for jet velocity 
ratios > 0.8. 

• Jets pointing directly upstream were ineffective due to weak vortex formation. 

• Jet arrays that cause counter-rotating vortex pairs can cause significant spanwise 
variations. 

• Spanwise average reattachment location seems to be more strongly affected than the 
detachment location 
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C. 2-D FAN PASSAGE SIMULATION 



The promising results of the UTRC experiment suggested the need for an evaluation 
of the concept in an actual transonic fan configuration. A full-scale simulation would be 
desirable, but such an experiment has the potential problem of not being definitive 
because of extreme difficulty in measurement and, therefore, in adequately assessing the 
results. Hence, a steady 2-D stream- surface simulation was undertaken in the present 
study. The combination of an experimental and computational simulation was the 
approach taken in an attempt to establish a representative baseline flow geometry in 
which to verify the effectiveness of the proposed vortex generators in alleviating the 
shock-induced separation. 

The geometry for the 2-D experiment was a simulation of the relative flow on an 
advanced fan rotor at approximately 63% of the span. The 2-D model used in the 
experiment was based on the stream surface geometry, but the blade profile was 
approximated (closely) as a wedge arc for ease in manufacture. This was a logical step 
since steamline contraction could not be simulated in the experiment. The geometry of 
the 2-D simulation is shown in Fig. 5. 

In reporting the generation of the baseline flow geometry in the present document. 
Section II describes the experimental simulation and Section HI describes the numerical 
simulation. As originally conceived, from the experimental simulation, measurements of 
surface pressures and schlieren photographs of the shock-boundary layer interaction were 
to be obtained. Measured boundary pressures were then to be used as input boundary 
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conditions for the numerical simulation, and a comparison could be made of both shock- 
boundary layer interaction structure and distribution of surface pressure. However, as 
described in the results given in Section IV, while the experimental apparatus was 
relatively successful, it could not be operated at sufficiently high back-pressures (>2x 
upstream) to produce a proper simulation of the fan shock structure at design conditions. 
Also given in Section IV are the results of a computation of the flow for the design 
pressure ratio of 2.28 using an Euler code. It was noted that, in the absence of viscosity 
in the code calculations at design pressure ratio, and with too low a back-pressure in the 
experiment, both computation and schlieren results showed qualitative agreement for 
shock structure at the entrance to the passage. It was concluded, in Section V, that 
improvement (through a refinement in the design of the valve) in the control of the back- 
pressure is required to properly simulate fan conditions. Grid refinement is required to 
continue the computational simulation. 

Details of the experiment are given in Appendices A - C and details of the 
numerical approach and the grid generation code are given in Appendices D - F. 
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II. EXPERIMENTAL SIMULATION 



A. TRANSONIC CASCADE WIND TUNNEL 

1. Wind Tunnel Description 

The original specifications, design and construction of the Transonic Cascade Wind 
Tunnel are described in Ref. 5. Subsequent modifications to the tunnel are described in 
Ref. 6. 

The apparatus is a blowdown device with a design Mach number of 1.4 and a 
design total pressure of 50 psia. A schematic of the Transonic Cascade Wind Tunnel is 
shown in Fig. 6. The tunnel was relocated (to Bldg. 216 at NPS) and operated with a 
new 300 psia air supply system in the course of the present work. A photograph of the 
right hand side of the wind tunnel is shown in Fig. 7, and a photograph of the left hand 
side is shown in Fig. 8. A tubular flow straightener is used to remove the swirl induced 
by the upstream pipe bends and the control valve. The settling chamber is composed of 
a flat, perforated plate diffuser followed by a screen to reduce the turbulence level and 
promote flow uniformity. There is a round (10 in. diameter) to rectangular (2.0 in x 4.5 
in) transition section from the settling chamber to the supersonic nozzle. The tunnel uses 
a 2-dimensional converging-diverging nozzle to accelerate to M = 1.4. A photograph of 
the interior of the test section (without blades) is shown in Fig. 9. Aft of the test section 
is a back pressure valve that is used to adjust the shock location within the test model. 
The back pressure valve is a "ramp and drum" throttle valve with the ramp actuated after 
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the flow "starts" by a pneumatic cylinder. A schematic of the back pressure valve is 
shown in Fig. 10, and a photograph of the exterior of the valve is shown in Fig. 1 1. The 
exhaust is ducted to a sound suppressor outside of the building. The tunnel start-up is 
controlled by a pneumatic control valve. A photograph of the control valve is shown in 
Fig. 12, and a photograph of the control panel is shown in Fig. 13. 

2. Data Acquisition System 

The data acquisition system used in the experiment is described in detail in Ref. 7. 
A schematic of the system is shown in Fig. 14, and a view of the system is shown in Fig. 
15. The BASIC code used to record the tunnel reference pressure measurements (using 
an uninstrumented bottom blade) is given in Appendix B, Fig. Bl. The Scanivalve port 
number and corresponding pressure location were as follows: 

• Port 1 - Atmospheric pressure 

• Port 2 - Calibration pressure (50 psia) 

• Port 4 - Plenum pressure 

• Port 6 - Upstream static pressure 

• Port 8 - Downstream static pressure 

Intervening ports were vented to the atmosphere to verify that the measured pressures had 
time to become steady state. 



9 



3. Optical System 

A schematic of the optical system is shown in Fig. 16. The system included a 
continuous light source schlieren system and a spark source schlieren system. The 
continuous light source was a 100 W mercury vapor arc lamp, and the spark source was 
a 200 - 300 nsec spark discharge. The parabolic lens provided a parallel light beam 
through the test section. The parabolic mirror, mirror 2 in Fig. 16, on the far side of the 
test section focused the parallel beam onto the knife edge located in front of the camera 
box unit. At the rear of the camera box was a viewing screen behind a shutter capable 
of automatic or manual operation. The minimum exposure time was one-thousandth of 
a second. A Polaroid film holder could be inserted into the camera box when needed. 
In the present experiment black-and-white Polaroid Type 52 film was used. A front 
surface flat mirror, mirror 1 in Fig. 10, was inserted between the continuous light source 
and the parabolic lens when the spark source was to be operated. Mirror one was 
positioned and aligned using a highly focused flashlight Both the continuous light source 
and the spark source were filtered to provide high resolution schlieren photographs. A 
view of the left hand side of the optical system is shown in Fig. 17. 

B. TEST SECTION DESIGN 
1. Design of the Blading 

Based on stream-surface conditions through the 3-D fan blade, the 2-D experiment 
was designed to retain the approximate geometry while simplifying curvatures for ease 
of manufacture. It was found that the stream-surface fan blade profile could be 
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approximated by a wedge-arc blade shape having a flat pressure side and equal leading 
and trailing edge radii as listed in Fig. 5. The leading and trailing edge radii were 0.015 
in, and the chord of the blade was 6.000 in. The suction surface starting at the leading 
edge was inclined to the pressure surface at 3.5° to a point 2.85 in. along the pressure 
surface. The suction surface wedge was blended to the trailing edge radius by an arc of 
radius 13.53 in. 

The test section geometry (oriented to allow the flow direction to be horizontal) was 
composed of three partial blades defining two blade passages, as shown in Fig 18. The 
top passage was to be the control passage while the bottom passage was to be used to 
study the effects of the vortex generators. For the experiment with the vortex generator 
jets, the bottom blade was for control, while holes were to be drilled in the center blade. 
A photograph of the model blade is shown in Fig. 19, while the machine drawing of the 
model blade is given in Appendix A, Fig. A 10. Figures 20 and 21 show views of the 
upper and lower partial blades respectively, while the machine drawings of the blades are 
given in Appendix A, Figs. A9 and All respectively. Appendix A, Fig. A12, details the 
location of the pressure taps for the instrumented lower blade. Figure 22 shows all three 
blades attached to the inner plate. 

2. Design of the Side Pieces, Inner Plates and Windows 

The design of the test section was driven by three major considerations: nozzle wall 
boundary layer diversion from the nozzle, wave cancellation from the leading edges of 
the blades and optical accessibility. Boundary layer diversion for the top and bottom of 
the test section was accomplished through bleed channels above the top blade and below 
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the bottom blade, using existing ducting [Ref. 5J. Side wall boundary layer diversion was 
provided by designing an inner plate which began near the leading edges of the blades. 
Figure 19 shows a photograph of the inner plate and center blade, while Appendix A, Fig. 
A2, shows the machine drawing of the inner plate. The exhaust for the diverted side wall 
boundary layer was through a slot in the side piece. Figure 23 shows a photograph of the 
side piece and Appendix A, Fig. Al, shows the machine drawing of the side piece. 
Appendix A, Fig. A4, details the design of the side piece exhaust slot. In addition, the 
side piece and inner plate incorporated a window for optical measurements. The window 
was made of 3/4 in. Plexiglas, and was designed to allow viewing access to the 
anticipated shock-boundary interaction regions in both passages. Preliminary tests 
indicated that the thick Plexiglas would give acceptable schlieren quality. Lastly, the 
cancellation of weak waves from the leading edge was effected by adjusting bleed through 
a porous wall on the top of the test section forward of the top blade [Ref. 6 and 8]. 

C. EXPERIMENTAL RESULTS 

Six test runs were made. The first experiment, Run 001, was with the back pressure 
valve initially completely open. The flow angle of incidence (to the suction surface) was 
set at 1.15°. Figure 24 shows the continuous light schlieren photograph which was 
obtained. The photograph shows an oblique shock leading from the leading edge of the 
middle blade to the suction surface of the bottom blade. The multiple waves on the top 
blade are from the porous wall, which was not adjusted to reduce reflections. The 
parallel weaker oblique shock in the lower passage is believed to originate at the 
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intersection of the leading edge of the supporting tab on the middle blade and the inner 
frame. Once the photograph was taken the back pressure valve was slowly closed in an 
effort to bring the oblique shock forward. Figure 25 shows the continuous light schlieren 
photograph when the back pressure valve was rotated to fully close the gap between the 
eccentrically mounted cylinder and the (fully raised) ramp. The oblique shock is seen to 
have moved forward towards the leading edge in both passages. The shock appears to 
be near normal impinging on the boundary layer. Tables IA and IB lists the pressure 
recordings for Run 001. 

The experimental pressure ratio p 2 /p, was found to be 0.765 with the back pressure 
valve open. With the back pressure valve fully closed, p 2 /p, = 1.39. Calculations for 
the fan geometry indicated that the pressure ratio should be 2.28 on design when a strong 
normal shock would exist across the blade passage. On inspection, the throat of the valve 
in the fully closed position when added to the spaces around the sides of the ramp gave 
an area which was too large to achieve the desired throttling. 

Run 002 was made after the actuator was adjusted to allow the ramp and drum to 
fully close when the drum was actuated. However, the pressure ratio p 2 /p, obtained in 
Run 002 was 1.26. No photographs were taken of this run. Run 003 was made with the 
actuation pressure increased to 600 psia. Again, inadequate pressure ratios were obtained. 
The schlieren photographs and pressure measurements (not shown) were very similar to 
those recorded in Run 001. 

Run 004, with the ramp actuator operated at 750 psia, showed no improvement. 
Run 005 was run with the back-pressure valve mechanically fixed in a "closed" position 
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using a pipe sleeve around the actuator rod. Of concern was the possibility of the flow 
not starting. Although the flow was started successfully through the model, it was found 
that the ramp was 1/16 - 1/8 in. lower than was required to fully close, and, again, the 
required control was inadequate. For Run 006, the ramp was wedged in the fully up 
position and the spaces between the ramp and duct walls were sealed using composite 
material plates. During this test, a maximum back-pressure ratio p 2 /p, = 1.98 was 
obtained before the wedge holding the pipe in place and one composite plate gave way. 
As the back-pressure approached 1.98, the shock location moved toward the leading edge 
of the blade; however, the motion was very unsteady and very sensitive to the back- 
pressure valve position. Table IC lists the pressures recorded in Run 006. No schlieren 
photograph was attempted in Run 006 in view of the need to observe the shock movement 
while experimenting with the valve position. 
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III. NUMERICAL SIMULATION 



A. COMPUTATIONAL SCHEME 

The computational code used in the present work was based on the thin-layer 
Navier-Stokes equations. The scheme used a body-fitted coordinate system, and was a 
full upwind algorithm based on Roe flux-difference splitting. The Baldwin-Lomax 
turbulence model was used except near shock impingement. The Johnson-King 
turbulence model was used to model the shock-induced separation. Lastly, the Osher- 
Chakravarthy TVD scheme for flux-limiting was used. Details of the development of the 
numerical code and background theory are found in Appendix D. 

B. GRID GENERATION 

The grid for the numerical scheme was generated using the GRAPE grid generation 
code, described in Appendix E. Two different grids were generated: a channel passage 
C-grid, to be used to carry out a preliminary calculation of the flow over a single blade 
between parallel walls, and a cascade C-grid, to be used to compute the flow through the 
simulated fan passage. 

The channel passage C-grid was an 199 x 64 grid suitable for Euler calculations. 
The input file for use with GRAPE is shown in Appendix F, Fig. FI. Figure 26 shows 
the channel passage C-grid. 
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The cascade C-grid was generated using a modification made to GRAPE which is 
outlined in Appendix E. The input file for GRAPE is shown in Appendix F, Fig. F2. 
The grid size was 169 x 40. The cascade feature was used to obtain the necessary flow 
periodicity in the developed grid. Figure 27 shows the complete cascade C-grid, while 
Fig. 28 shows an enlarged view of the leading edge region. Figures 29 and 30 show the 
details near the leading and trailing edges respectively. Figure 31 shows two cascade C- 
grids plotted simultaneously with one grid at an offset equal to the blade spacing. This 
figure demonstrates the periodicity of the grid and the accuracy of the geometric 
specifications such as the solidity. Figure 32 shows a detailed view of the boundary 
between the two grids. 

C. COMPUTED RESULTS 

The Mach number contours for the Euler solution obtained for the flow over the 
single blade in a channel using the channel passage C-grid are shown in Fig. 33. The 
solution shown in Fig. 33 is not a steady state condition but rather a "snapshot" of the 
flow during a transient in which the pressure side shock was propagating towards the 
leading edge. 

The Mach number contours for the Euler solution for the flow through the cascade 
C-grid are shown in Fig. 34. The angle of incidence of the flow to the suction surface 
was 0° and pj/p, = 2.28. The contours show a pressure side leading edge oblique shock 
and a reflection from the suction surface, followed by a strong, near-normal passage 
shock. Figure 35 shows the pressure contours. On both sets of contours, the leading 
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edge disturbance does not continue to propagate smoothly across grid boundaries because 
of a discontinuity encountered there in the grid spacing (see Fig. 32). The calculated 
blade surface C p distribution is shown in Fig. 36. This figure also demonstrates the 
effectiveness of flux-limiters in suppressing the oscillations that develop around the 
discontinuities. 
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IV. DISCUSSION OF RESULTS 



A. EXPERIMENTAL 

The oblique shock patterns within the passage, shown in Figs. 24 and 25, were 
observed to be qualitatively similar to the computed shock pattern shown in Fig. 34. It 
is noted, however, that the flow incidence to the suction surface in the experiment was 
1.15° whereas it was 0° in the computation. The shock should therefore have been 
stronger in the experiment than in the computation, and this certainly was the case. 

In the present work, the difficulties with the back-pressure valve precluded obtaining 
a pressure ratio p 2 /p, = 2.28. However, considerable experience was gained in the 
attempts which were made to achieve this ratio. First, it was found that the flow could 
be started through the test section with the ramp fully up. This means that the ramp can 
be fixed rigidly to the duct wall and need not be supported by a pneumatic actuator. The 
rigid positioning will allow the sides of the ramp to be sealed completely and will 
eliminate a potential source of unsteadiness in the valve throat area. 

For the brief period during Run 006 that the back-pressure was held at a value 
approaching twice the inlet pressure, the flow was highly unsteady. Unfortunately, the 
inlet pressure was dropping as the supply pressure fell, and it was not possible to 
experiment further. It was noted on the schlieren screen that the lower passage unstarted 
while the upper passage remained started, raising the possibility that the achievement of 
similar flows in the two passages may not be automatic. 
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B. NUMERICAL 



The Euler solution using the cascade C-grid gave useful qualitative information on 
the structure of the flow field, which was supported by the experimental data. It is 
expected that at the same pressure ratio, the viscous solution should move the shock in 
the passage due to losses and the displacement effects of the boundary layer. The flow 
structure might still not correspond to what is expected to occur in the fan, wherein the 
normal shock sits at the leading edge in an "unstarted" structure. It is also possible that 
the computed structure is the result of the initial conditions imposed on the unsteady 
solution. What is found is a solution that requires a "starting" process that is not present 
physically in the fan application. What would be required to obtain the "unstarted" 
solution would be to first increase the pressure ratio to expel the shock, and then decrease 
the pressure ratio to the design value. It is interesting to note that this is exactly the 
procedure which must be followed in the experiment, which is why the back pressure 
control valve is such a critical component in the test apparatus. 

Grid refinement is also necessary for the intended application with a viscous flow 
solver. The boundary layer will need to be resolved adequately with an appropriate 
number of grid points. Also, the orthogonality between the upper and lower boundaries 
of the cascade will have to be enforced more rigidly (Fig. 32). Lastly, the grid spacing 
will need to be more uniform to ensure continuity across boundaries (Fig. 32). Continuity 
in grid spacing is expected to yield more accurate propagation of shock discontinuities. 
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It is noted, however, that the validity of the turbulence model chosen has been 
demonstrated (Appendix D, Fig. Dl), and a solution for an airfoil has been obtained using 
the viscous solver (Appendix D, Fig. D2). 
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V. CONCLUSIONS AND RECOMMENDATIONS 



In the present study, both experimental and numerical simulations of the flow on 
a stream surface through a transonic fan passage were attempted with a view to providing 
a baseline geometry in which the effectiveness of upstream vortex generating devices in 
alleviating effects of shock-boundary layer interaction could be examined. From the 
experiment, the following were concluded: 



• The design of the test section geometry was successful. No structural problems 
were encountered in six test runs. Boundary layer diversion slots appeared to start 
by natural aspiration to the atmosphere. Schlieren photographs were obtained 
successfully through thick Plexiglas windows. 

• The upstream pressure control valve was fully satisfactory, and the transonic 
cascade wind tunnel was made operative in its new location. 

• The downstream pressure control valve was found not to be satisfactory as 
originally installed and was shown to be critical to the eventual success of the 
experiment. It was shown that the test section flow would establish with the valve 
ramp fully up, so that the ramp could be made both rigid and sealed. The eccentric 
aluminum cylinder component did not provide adequate control because of internal 
leakage. 

• No problem other than the back-pressure control valve was found which would 
prevent a successful experiment. 

From the initial steps taken toward a numerical simulation, the following were 
concluded: 
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• A suitable grid for computing the cascade flow using an Euler solver was generated 
using GRAPE code as modified by R. Chima (Appendix E). A suitable grid for 
viscous calculations near this blade surface can be obtained using this code. 
Discontinuities in grid spacing across periodic boundaries given by the code can 
possibly be removed using the DSOBI parameter (Appendix E). 

• The solution obtained for the experimental cascade at the design pressure ratio 
using an Euler solver showed a normal shock wave well inside the passage. Since 
the inclusion of viscous displacement effects should move the shock downstream, 
this solution may be the result of the assumed initial conditions. 

• The oblique shock structure predicted within the passage was qualitatively similar 
to that observed in the schlieren photographs obtained from the experiment at lower 
pressure ratios. 



The following are recommended to advance the experiment: 



• Perfect the back pressure control valve by sealing and rigidly attaching the ramp 
and replacing or redesigning the eccentric drum. 

• Install pressure taps in the -lower blade (as designed) and connect to a high-speed 
scanning data system. 

• Manufacture additional windows, or provide aluminum blanks for tests not requiring 
optical access. 

• Provide a tight shut-off valve in series with the present tunnel control valve. 

• Provide a safety rupture diaphragm upstream of the back pressure control valve. 



The following are recommended to advance the numerical simulation: 



• Ensure grid spacing continuity, orthogonality along the boundary of the grid and 
improved resolution in the boundary layer region for use with the viscous solver. 

• Investigate the "starting" process. 

• Run the Euler solver with the flow incidence at 1.15° to the suction surface. 
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TABLE IA. Run 001 Pressure Data, p 2 /Pi = 0.765 



Port 


Pressure (psia) 


1 

2 

4 

6 

8 


14.8 
50.0 
40.3 

17.9 
13.7 



TABLE IB. Run 001 Pressure Data, 


P 2 /P 1 = 1.39 


Port 


Pressure (psia) 


1 

2 

4 

6 

8 


14 . 8 

50.2 
22 . 4 

14.3 

19.8 



TABLE IC. Run 006 Pressure Data, 


P?/Pi = 1.97 


Port 


Pressure 


1 

2 

4 

6 

8 


55.8 

18.1 

35.7 
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Figure 1. Shock-Induced Separation on a Fan Blade [Ref. 1] 
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Figure 3. Typical Surface Pressure Distribution 
in the Interaction Region 
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Blade Geometry 



L.E. Radius = 0.015 in 

T.E. Radius = 0.015 in 

Wedge Angle = 3.5° 




Figure 5. Cascade Geometry 
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Figure 6. Schematic of Transonic Cascade Wind Tunnel 
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Figure 7. RHS of Transonic Cascade Wind Tunnel 
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Figure 8. LHS of Transonic Cascade Wind Tunnel 
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Figure 9. Test Section Interior (w/o blades) 




Figure 10. Schematic of Back Pressure Valve 
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Figure 11. Back Pressure Valve 
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Figure 12. Control Valve 
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1. HG-78K Scanivalve Controller 

2. HP 3455A Digital Voltmeter 

3. HP 3495A Scanner 

4. HP 9862A Calculator Plotter 

5. HP 9830A Computer 

6. HP 9868 I/O Expander 

7. HP 11305A Controller 

8. HP 9867B Mass Memory 

9. HP 13215 Disc Power Supply 

10. Pressure Transducer (48 Port Scanivalve) 

11. HP-IB Interface Bus 



Figure 14. Schematic of the Data Acquisition System 
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Figure 15. Data Acquisition System 
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Figure 16. Schematic of Optical System 
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Figure 17 



LHS of Optical System 
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Figure 18. Schematic of Test Section 
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Figure 19. Test Section Inner Wall and Center Blade 
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Figure 20. Test Section Upper Blade 
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Figure 21. 




Test Section Lower Blade 
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Figure 22. Test Section Inner Plate and Model Blade 
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Figure 23. Test Section Side Piece 
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Figure 24. Run 001 Schlieren - P 2 /P 1 - 0.765 




Figure 25. Run 001 Schlieren - p 2 /pi = 1.39 
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Figure 27. Cascade C-grid 
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Figure 36. Cascade Blade Surface C p Distribution 
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APPENDIX A 



MACHINE DRAWINGS OF TEST SECTION COMPONENTS 



Drawing No. 


Title 


T001 


Test Section Side Piece 


T002 


Test Section Inner Frame 


T003 


Test Section Inner Frame and Blades 


T004 


Detailed View of Side Piece Slot 


T005 


Side Piece Window Retainer Frame 


T006 


Side Piece Window 


T007 


Sectioned View of Window Installation 


T008 


Detailed View of Window Slot 


T009 


Test Section Top Blade 


TO 10 


Test Section Middle Blade 


T011 


Test Section Lower Blade 


T012 


Pressure Tap Locations on Lower Blade 
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Figure A1 . Test Section Side Piece 
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Figure A2 . Test Section Inner Frame 
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Figure A3. Test Section Inner Frame and Blades 
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Figure A5. Side Piece Window Retainer Frame 
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Figure A6. Side Piece Window 
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Figure A7 . Sectioned View of Window Section 



64 



o ° 
o 

4 l 




Figure A8 . Detailed View of Window Slot 
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Figure A9 . Test Section Top Blade 
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Figure A10. Test Section Middle Blade 
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Figure All. Test Section Lower Blade 
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APPENDIX B 



DATA ACQUISITION PROGRAM 



10 REM *-***17:C001***** F'.P.SHF'EEVE 14 *1 

20 REM 

20 REM THIS P POGROM IS FOP DATA ACQUIS I TInN 
TO REM FROM A SCflMI VALVE ON THE TRANSOM If rASrRDE 
50 REM 

60 REM VARIABLES: 

70 REM 
80 REM 
90 REM 
100 REM 
110 PEM 
120 REM 
130 REM 
140 REM 
150 REM 
160 REM 

170 REM DATE: MAR 14 < 1 9?A 

ISO REM 

ISO DIM RE 5* 10 3* DC 5* 10 3* PC 10 3* P*[ 20 3* A-TC k J 
200 MAT A=ZER 
210 MAT D = ZEF: 

220 MAT R=ZER 

230 PE 1 3=1 

231 PC 2 3 = 2 

232 PE 3 3 = 3 
234 P[ 4 3=4 

• 235 RE 5 3=5 

, 237 RE 6 3=6 

233 F E 7 3=7 
23? PE 8 3=8 

! 240 DISP "EATER ATMOSPHERIC RRES C = Mil HG • " 5 

i 250 INPUT F'8 

j 260 DISP "ENTER NO. POINTS TO TAKE (MA^=5- M : 

! ’270 INPUT SI 

275 DISP "PRESS <CONT> WHEN READY" 5 
| 276 STOP 

I 2S0 FORMAT B 
290 FORMAT 2B 
: 300 FORMAT 4B 

! 310 FORMAT F3.0 

320 REM 

330 REM INITIALIZE BUS AND ASSIGN DATA ELEMENT?; 

340 REM 

350 WRITE < 1 3 !* 300)256* 20* 768*512 5 
360 CMD "?D#"p "F1R7M3R0H0T3" 

370 REM 
330 PEM 
390 REM 
400 V= 1 
410 C 1 =0 
420 S= 1 
430 REM 

440 REM BEGIN DATA SAMPLING 
450 REM 

460 FOR K=1 TO SI 
470 FOR J=1 TO 8 



Figure Bl. Data Acquisition Program TXC001 



70 



4 77 


am 


D 4 




4 76 


WRITE 


' 1 3 - 2 


80 > V «' 


4 80 


A = FT J] 






4 80 


GUSUB 


1470 




700 


T 1=0 






710 


PPM 






720 


Pin PF'ESSUP 


E LINT A 


7 30 


PPM 






740 


am 


IN " 




770 


WRITE 


< 1 8 * 2 


80 ’ V 5 


760 


FAR L = 


1 TO 


er 


770 


am M ? 


D ! ” 




780 


Wf'IIE 


•13*3 


1 0 ' V + 8 


780 


am M ? 


LOT * " 


T 3" 


600 


cud “••• 


l r " 




6 1 0 


ENTER 


as* * 


•Dl 


620 


T 1 = T 1 + 


Dl 




6 30 


me::t l 






6 4 0 


Ml =T 1 


er 




670 


dl r * j : 


i=ni 





67 1 JH}*10]=P8 

660 mf::t j 
f. 7 n r 5 i 1 1 n 77n 

PEN P E M Li AINUS TEMREPATUFE OH SCANNER CHANNEL 0 

680 REM 

700 Ctm “ ? Li < ” * " C “ 

710 Ml.ll PUT a 3- 310 'Ll 

720 am ■"■•Dtr - m t. 3 m 
730 am M ?ci" 

740 ENTER • 13. *H'C t . 7 ] 

770 C nil M ?IK M « "C" 

760 PEN 
770 PEN 
780 PRINT 

780 PRINT " PUT A SAMPLES FOP POINT " f 

800 PRINT 

810 FOP 1=1 TO 10 

820 WRITE a 7*830.’ I .HU * I T 

8 30 rORMHT F3. 0.7::* F 13.6 

840 me::t i 

870 rasp "IS DOT ft VAL I D? 1 = VES 0 =F:EThLE 

860 INPUT h2 

870 IF Fl2= 1 THEN 880 

880 GOTO 470 

880 next k 

800 PEN 

810 REM DEFINITION OF DOTH ftRPftV 
820 REN 

840 REN LK*?1> = OPEN •; ATMOS PRESS. > 

870 REN IK*'2) = CALIBRATION PRESS. (SUIT CM ED > 

860 REN D a * 3 > = OPEN 

870 PEN IK * ? 4 > = PLENUM PRESSURE 

880 PEN IK*. 70 = OPEN 

880 REN DC*? 6 > = UPSTREAM STATIC PRESSURE 
1 000 REN D < * ? 7 > = OPEN 

1010 PEN IK * ? 8 > = DOWNSTREAM STATIC PRESSURE 
1020 PEN DC** 8 > = OPEN 

10 30 PEN DC*? 10 > = ATMOSPHERIC PRESSURE 
1070 PEN 

1060 REN DAT ft CONVERSION l STORAGE 
1070 REN 

1080 FOR » = 1 TO SI 
1080 AIK* 1 ]=DU:? 1 ]- 100000 
1 100 At K* £ ]=DC \ * £ ]* 100000 
1110 AC K ? 3 ] = DC V ? 3 1 00000 
1 120 ACK*4 ] = Dt K ? 4 ]* 100000 



Figure B1 (cont) . Data Acquisition Program TXC001 
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1 1 30 fin • 5 3^- lit 1 . 53*1 00000 
1 1 40 I’ll I . • 6 J~li[ I . 6 3- 1 00*300 

1100 nn:. r 3- in i > r 3 * 1*300*34. 6i4=> 

1 1 6 0 I’ll 1-8 3 — I* 1 1 . 8 3 

1 iro m 1 • ■? 3 ~ 11 1 1 . ? .1 

I 18*3 I'll I . 1 1 . 1 3 = 1*1 I • 10 3 

II oo im:i 1 

11 -?5 11 1 si • " open iirn n file nun 1 : 11111 . 12 * 10 "; 

1196 81 Of 

1200 1USP "Eli I ER I'll r H TILE III HIE"! 

1210 mrui nr 
1220 tiles * 

12 30 038 1 CM flt.l.l? 

1240 Miir n un tt m* 

1 250 r F : I II 1 "HfiTfi 81 OPED 111 THE "fit 
1260 PR I 111 

12 re F'EM 

14 30 PRIlir 

1436 PRINT "fll HOSPIIEP I C PRE08I IFF = "p:3 
1440 Ellli 
1*1 5 U 

1460 F'EM SUBRuUt I HE "POSIT" 

147*3 GOSUB 165*3 
1480 H=n- P 
14? 0 CUD " ’ li!" 

1500 IF H 0 1IIEII 1580 
1510 IP H o THEM 158*3 
1520 RETURN 

1530 FEM home s v 

1540 I IF 1 1 E <: 1 3. :31*3 .y +4 
1550 HR I IE ( 1 3 • » ' “ P " 

1560 Ill’ll I 4*30 

157*3 60 fo 147*3 

158*3 REM flUVflHCE 3 V 

1 6 *3 *3 I IF: HE <13.31 *3 > V - 1 

1610 HR I IE < 13.* ••T" 

1620 unn 50 

1640 GOTO 1470 

1650 FEM PERU 8 V FWDPESS 

1660 CMD "?GT " 

1670 F'0=RBYTE 1 3 

1680 l=ei i null (pons:' 

1690 T=F:OT-:R0.4) 

1700 M=P I nilLi • T . 7 :■ 

1710 F'= 1 0*M + L 

1720 HR HE <13. 28*3>256. 8=;; 

1730 RETURN 



Figure B1 (cont) . Data Acquisition Program TXC001 
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APPENDIX C 



PROGRAM TO DRAW A WEDGE-ARC CASCADE 



10 PEN *****c;3CB03*-****- F’ . F . SURE EVE 

£0 PEN PROGRAM TO CflLC. ? DR Fill R HEDGE-ARC CASCADE 

30 PEN INTO n cascade 

35 DEG 

40 uniR 6.Q.015«3.5 i2.85*-51.84i2*?0 
50 PERU CO » RO » A 1 « M 1 » GO • SO « G? 

54 A0=-G0 

55 D 1 = P 0 * r 1 -C OSAO-8 I NmO > 

5S D2 = F'0« I-COShO+SIHAO) 

57 S1=C0*SINAQ 30 

53 S2=CO SO 

SO PEN CRLC . BLADE GEOH. 

70 «o=PO 



30 VO=RO 



1 00 


j;i =xo- 


po*smni 


i io 


Yl=PO*R0«tOSfll 


120 


N2-R0 




130 


‘.'2=0 




140 


E0 = < F:0 


i hin hi • 2; ;-::2 


150 


A8=C0- 


PO-Hl 


ISO 


R9=A8* 


HP 


170 


B8=F:0- 


' H1+E0>*TAIIR1 


180 


B9=B8* 


38 


130 


p l = ( R9+B9-R0*R0 > / < 2 * * R3 ! 


200 


;:?=ni 




210 


V3*<Hl+E0)*Tfllini 


220 


?M=X '3 + R 1 ^smni 


230 


Y4=Y3- 


■pi*cqsai 


240 


:'.7=co- 


F:0 


250 


V7 = F'0 




2 SO 


r i =mt 


K <5<7-JM >'< Y7-Y4 > > 


2 70 


'{5=X r + F 0*3 IMP 1 


230 


V5=V7+RO*COSP 1 


2?0 


X6=X7 




300 


'v k - n 




310 


PEN — 


PR IIIT RESULTS 


320 


PRINT 




330 


PRINT 




340 


PRINT 


” HEDGE -ARC CfiCHDE 


350 


PRINT 




3S0 


PRINT 


"CHOPD="CO« "HEDGE 


370 


PRINT 


"HEDGE HNGLE= " H 1 " 


380 


PRINT 


"SURFACE RADIUS*" 


3 90 


PRINT 




400 


PRINT 


"xi= M .'{i» •• v i = •• v i 


410 


PRINT 


" J{2="}{2j " Y2=" V2 



-PO ■- 



BLADE " 

LEIIGTN="U1 

DEG *' » " L . E . F'MD 1 03 = *' PO 
R 1 » *’ F'H I = ” P 1 " DEG. " 



10 1? ?u 
BLADE 



Figure Cl. Cascade Geometry Program SSCB03 
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420 
4 30 

440 

450 
4 £ 0 

470 

471 

472 
4 75 
4 3 0 
4 90 
1 000 
i o i 6 
1 020 
1 020 

1040 

1041 

1042 
1 050 
1 0£0 
lObi 
1 0£ 2 
1070 
1 020 
1021 
1022 

I 090 
1100 
1110 
1 1 1 1 

I I 12 
1120 
1 1 20 
1 1 31 
1122 
1 140 

I 150 
1151 

I I 52 

1154 

1 1 55 

1156 
1160 
5000 
5010 
5020 
5025 
5030 



"yA = ‘ 
"X 5=‘ 
M XS=' 



M- 

■'5 - 

n 



V4 = ' 



• i :• 

V4 

Y5 

Y6 



•cnscniiE 

•set ting 



GEOMETRY: " 
HMGLE ” " 00 " L 1 E G . 



FOP SUL 1 1' I T Y = " 30 * *' BLOI'E SFmC I NG = " S2 



r F I M T 
rr-:iMT 
PRINT 
PRINT 
rpiMT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
STOP 

PPM DPiUlPLnDE 

SC OLE - 1 * 14? 2-3 
LllBEL <*-1.5-2-0-10 150 
TOP T = -90 TO 90-fii STEP 
X9 = R0 * < 1 -C US ] } 

Y9=P0*C 1+SIN1 ' 

GO SUB 5000 
NEXT T 



•BLNlif 3rHCE<N.' = ,, Sl 



130-01 



30 



Y9 = Y3 

GO SUB 5000 

FOR T=-M1 TO FI STEF <nifFlv250 

::?=:m+ri*sint 

Y9=Y4+R 1 * LOST 
GO SUB 5000 
NEXT T 

FUR T=-90*F1 TO 90 STEP < 180-PI >30 

X9=X7+R0*C0ST 

Y9=Y7-F0*SINT 

GOSUB 5000 

NEXT T 

X9=F:Q 

Y9=0 

GOSUB 5000 
FEN 

DISP " INPUT 2 FOP BLOUE 2* 3 FOR 3 " : 

INPUT ZO 

0 1 = < ZO - 1 > * 3 2 * C0SG9 
02= <20-1 >^S2*S I NG° 

OFFSET 01-02 
GOTO 1030 
ENH 

PEN SUBROUTINE TO TF0N3L0TE ONH ROTATE FIXES 

X = < X9 - II 1 > * C 0 S G 0 + < V 9 - U 2 > * S I NG 0 
Y= < Y ?- D2 > *C0SG0 - < X?- D 1 > *S I NGO 
PLOT X* Y 
RETURN 



Figure Cl (cont) . Cascade Geometry Program SSCB03 
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APPENDIX D 



DESCRIPTION OF THE NUMERICAL METHOD 
A. REYNOLDS-AVERAGED NAVIER-STOKES EQUATIONS 

The Reynold’s-averaged Navier-Stokes equations are presented first. The unsteady, 

compressible Navier-Stokes equations are developed in Ref. 9. A complete derivation of 

the Reynold’s averaging process can be found in Ref. 10, which is the source of the 

following description. 

The time-dependent, compressible Navier-Stokes equations in tensor notation are 
as follows: 

Pt + = 0 

(pu i ) t + (p UjUj + 0^)^ = 0 i = 1,2,3 (l) 

(p E) t + (p EUj + UjOjj + = 0 

where subscript notation has been used to denote partial derivatives, eg. p t = and the 

summation convention has been used for repeated indices. The stress tensor for a 
Newtonian fluid is given as 

°ij = + u j,l ~ j u k.k) < 2 > 

Here the Stokes hypothesis 3X + 2 =0 has been used. The heat flux vector is obtained 
from the Fourier law of heat conduction as 
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( 3 ) 



* - kT j = -p- h j 

X 

To derive the Reynold’ s-averaged Navier-Stokes equations, the averaging operation 
is defined as follows: 



t + T 

P (Xj,t)=jj,fp(Xj,S)dS (4) 

1 C-T 

where 2T is the time averaging interval, which is assumed to be large enough compared 
to the energy containing turbulent time scales but small when compared to the time scale 
of the average motion. It is assumed that there is a range of values of T where the mean 
density is independent of T. 

The fluctuating density is defined as the difference between the density and its mean 

value. 



= p - p (5) 

The averages and fluctuating values for the other variables are defined similarly. It is 
assumed that the average of the fluctuations equals zero and the average of the average 
still is the average. 

The time-averaged Navier-Stokes equations are obtained by averaging Eq. (1). The 
resulting equations involve terms similar to the following: 
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( 6 ) 



pu i u j/ - P U t Uj + p U' 1 u' j + U 1 P'Lf'j + Uj p' U' 1 + p ' u’ t U' J 

For incompressible flows, where the average of the density fluctuation equals zero, 
the last three terms of Eq. (5) are dropped. The compressible equations contain additional 
terms, such as it"’" p 1 u ~’ , therefore, an alternative method for averaging velocity and 
energy terms is followed. This method is called mass-averaging (Favre-averaging) with 
the following definitions: 

Qj = puf/p h = p7i/p < 7 > 



u' i = u i -u i h' = h - h (8) 

It is important to note that the averages of the fluctuating quantities for Favre-averaging 
is not zero. 

u ' ' j = -p ' u ' j/p h' ' - -p' h' /p ( 9 ) 

However, the mass-weighted variables are zero. 

pit' ' j = 0 p TT' ' = o (io) 

By introducing the mean and fluctuating quantities into Eq. (1) and averaging, the 
Reynold’s-averaged Navier-Stokes equations are obtained. 

p^ + (puj) j = 0 



(P uj t + (p + o^) j = 0 



(ID 
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(p E) t + (p EQj + UjOjj + qj) 



°ij = °i/ + °i/' 



The following are identified by name: 



Reynold's Stress Tensor: 



j = 0 i= 1,2,3 

<3j = Qj" + Qj 7 






r 



p Uj ' ' Uj ' ' 



o 

Viscous Stress Tensor-. a i:j M - + Uj :i - —t> 1 ;j u kk ) 

Reynold's Heat Flux Vector: q^ = p TP Hip ' (12) 

Molecular Heat Flux Vector: q, M = 

j Pr J 



Total Energy: E = e + k + u k u k /2 



Turbulent KE: k = pu*.' ' u k ’ ' /2p 

It is assumed that |i is independent of time, but depends on temperature variations. 

The Reynold’s stresses are related to the mean flow quantities through turbulence 
modeling. The turbulence models used in the present work are discussed in the next 
section. When the Reynold’s stresses and heat fluxes are related algebraically to the 
mean flow quantities, the resulting models are called algebraic turbulence models. The 
algebraic turbulence models are either mixing length-type models or eddy viscosity-type 
models. In the present study the widely used Baldwin-Lomax algebraic eddy-viscosity 
model was used. 
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The eddy viscosity models are zero-equation models, and the eddy viscosity is 
modeled algebraically. More complex models exist that take into account the transport 
of certain turbulence quantities such as kinetic energy or shear stress. These models are 
derived from simplified forms of the Navier-Stokes equations for mean and fluctuating 
quantities [Ref. 11]. 

A simple one-equation model that is specially designed to predict shock-boundary 
layer separated flows, the Johnson-King model [Ref. 12] was used in the present 
investigation. This model was considered to be most suitable for the computation of the 
nonequilibrium flow in the recovery region after the shock-induced separation, based on 
the experimental data given in Ref. 12. Figure D1 shows a comparison of predictions for 
the separated flow region after a transonic shock using four turbulence models. The 
Johnson-King and RNG models exhibit the best results. The Baldwin-Lomax model and 
the Johnson-King model are briefly described in the next section. 

B. TURBULENCE MODELING 

1. Baldwin-Lomax Model 

The Baldwin-Lomax model is a two layer algebraic model which does not require 
the finding of the boundary layer quantities, as does the Cebeci-Smith model. Instead, 
the outer and inner eddy viscosities are scaled by the vorticity. A complete formulation 
of the model is given in Ref. 13. 
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2. Johnson-King Model 

The shock- separated region is computed with the Johnson- King turbulence model. 
A detailed description of the development of the model is given in Ref 12. The basic 
formulation of the model necessary for the numerical implementation is presented briefly. 

The Johnson-King turbulence model was developed for the purpose of improving 
the predictions of pressure driven separated flows and transonic shock-induced separated 
flows. The Cebeci-Smith model [Ref. 10] was proven fairly accurate, even though it was 
based on the invalid assumption that the turbulent shear stress depends only on the local 
properties of the mean flow. The Johnson-King model develops a closure for boundary- 
layer flows without attempting to predict the production, dissipation and diffusion of the 
turbulent kinetic energy and the dissipation rate, throughout the entire viscous layer, as 
in two-equation models. 

To account for convection and diffusion effects on the Reynolds shear stress 
development, an ODE is developed to describe the maximum Reynolds shear stress 
development in the streamwise direction. The ODE describes the development of 
it v r = JxJp along the path of maximum shear stress. 

The eddy viscosity is given by 

v £ = v t [1 - exp ( — ) ] (13) 

where v u and describe the eddy viscosity variation in the inner and outer pans of the 
boundary layer. The inner eddy viscosity is computed as 
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v tj = D 2 kzJ(^u' v' ) MX (14) 

D 2 = 1 - e (_z/;r) 

where the constant A + = 15. The outer eddy viscosity is given by 

v t = a(x) -[0.0168C7 e 6*Y] < 15 > 

where y is Klebanoff s intermittency function 

y = [1 + 5 . 5 (4) 6 ] _1 (16) 

and o(x) is the solution of the ODE mentioned earlier. 

To complete the formulation, an equation is needed for i r v'„. This equation is 
developed from the turbulence kinetic energy equation. The result is the following ODE 
after a change of variables: 



$9 - r (2 9_\ + C dif L rn h _ ( Vt ° ) U] (17) 

dx 2^ lU gj ai 6[0.7 - (y/5)J 11 ( v > |J ‘ 



where Qlif and a, are modeling constants, u B is the maximum average mean velocity, g 
and g^ are given as 



g = 1 /y (-u* V )„ and = l/yj (-u* v* ) “ ~ (18) 

and L,,, is the dissipation length evaluated as 

L m = 0.40z zj5< 0.225 (19) 

1^ = 0.096 zjd > 0.225 

The equilibrium shear stress, g^, is determined from the following equilibrium eddy 
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viscosity distribution 



V t,e<7 V £„ 



eq 



[1 



exp (- 



V Cj . eg ^ j 
v 

c 0 , ecr 



ti.eq 



- D 2 kzJ(-u' V ) 



m, e<7 



( 20 ) 



V C 0 , e<y = 0 .0l68u e 6*Y 

An implicit Euler method is used for the numerical solution of Eq.(17), and the 
maximum shear stress at each iteration level is updated as follows: 



a (x ) n+1 



o (x) 




( 21 ) 



Solutions with the Johnson-King turbulence model were obtained as follows. First a 
convergent solution using the Baldwin-Lomax turbulence model for the entire flow field 
was obtained. Then the Johnson-King model was applied to the upper part of the airfoil. 
To initiate the solution o(x) in Eq.(15) is set equal to one, and it is allowed to change 
according to Eq.(21) until the final solution is obtained. Note that the Johnson-King 
model reduces to the Cebeci-Smith model when o(x) is identically equal to one. 



C. COMPUTATIONAL METHOD 

The following section was developed using notes from J. Ekaterinaris. 

1. Discontinuous Solutions and Entropy Method 

The numerical solution of the initial value problem for hyperbolic systems of 
conservation laws consists of the solution of 
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U t + f x {u) = 0 



( 22 ) 



u {x, 0 ) = Uq {x) 

where u(x,t) is the vector of m unknowns and f(u) is the flux, a vector valued function 
of m components. The matrix form of Eq.(22) is 

u t + Au x = 0 

Because the system is assumed hyperbolic, all the eigenvalues of the Jacobian matrix are 
real and positive. 

To allow for discontinuous solutions we admit weak solutions which satisfy Eq.(22) 
in the sense of distribution theory as shown below 

« M H 

J J {4> t u + <}) x f (u) } dxdt + J* 4> (x, 0) u 0 (x) dx = 0 < 23 > 

z- 0x=-«* -» 

where <)) are C*° test functions that vanish for |x| + t large. The hyperbolic system is 
satisfied in a pointwise sense by a piecewise smooth solution in each smooth region, 
while across each discontinuity the Rankine-Hugoniot relation 

f(u R ) - f{u L ) = C{u R - u L ) (24) 

holds where C is the speed of propagation of the discontinuity, and the left and right 
states are designated u L and u R , respectively. 

The conservation law Eq.(22) is considered to possess an entropy function E(u) 
defined as 
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• E is a convex function of u, E uu > 0 



convexity 



• E satisfies E u f u = F u compatability 

where F is the entropy flux function. 

The convexity condition forces irreversible processes to run in the correct direction 
and produce entropy. A consequence of the compatability condition is that the reversible 
processes do not produce entropy. Note that the mathematical entropy E, and physical 

i 

entropy S = log(p/p T ), statements are compatable, and if one is true then the other is true. 
Every smooth solution of Eq.(22) then satisfies 

E U (C) + F x (u) =0 (25) 

If u is piecewise smooth with discontinuities, then Eq.(25) holds pointwise in the smooth 
regions; and across a discontinuity 

F(u r ) - F(u l ) z C[E(u r ) - E(u l )] (26) 

Consider for simplicity that the numerical approximation to weak solutions of 
Eq.(22) are obtained by the following explicit scheme 

u i 1 = u i ~ "f” (f 1 * 1/2 ~ f 1 - 1 / 2 ) (27) 

where fu 1/2 is a numerical flux defined as 

f 1 . 1/2 = f(u?,u? +1 ) 

It is required that the numerical flux is consistent with the physical flux in the following 
sense 
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f{u, u) = f(u) 



In addition, the differencing scheme is consistent with the entropy condition if the 
inequality of 

£? +1 * Ei - + Fi- 1/2) < 28 > 

is satisfied where 

E I* 1 = EiuT 1 ) 

and 

Fi. 1/2 = , Ui n +1 ) 

Assume that the difference scheme in Eq.(27) is consistent with the hyperbolic 
conservation law Eq.(22) and the entropy condition Eq.(26). Let u° be a solution of 
Eq.(27) for initial values u° = O(j'A) . Suppose that for lim Ax — > 0 (Ax = x/i^) and At 
-> 0 (At = T/N) 

limAx-o uf = u{x,t) 

At-0 

the limit exists in the sense of L, norm convergence. Then the limit satisfies the weak 
form of the conservation law Eq.(23) and the weak form of the entropy condition given 
by 
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-/ J (4> C E + <t>j <E) dxdC - J 4> ( x , 0 ) E(u 0 (x) ) dx £ 0 
0 -«• -«• 

where <}) is a nonnegative smooth test function <t>(x,t) possessing a compact support. 
Assume that the initial data corresponds to some reference state u’ for large |x| 

u 0 {x) = u* |x| > M 



then 



Ui = u* for A xi > M + nAx 

Without altering the convexity of the entropy function an inhomogeneous function can 
be added to obtain 

E(u m ) = 0 E Ui (u m ) = 0 
By summing Eq.(28) with respect to i to obtain 

E 1 * E B " 

1 i 

shows that entropy is a decreasing function of time, and 

E E ° * E E ° 

l 1 

indicates that the scheme is stable. 

This condition is not strong enough to prove pointwise boundness of the solutions 
or the existence of convergent subsequences. Solutions of equations with compressible 
flow must ensure that the difference scheme used keeps the variables within their physical 
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range (density and pressure are always positive quantities). The above development is 
true for any number of space variables. In addition, they hold for multidimensional 
schemes consisting of one-dimensional factional steps provided that each individual one- 
dimensional step satisfies the entropy inequality. 

2. Upwinding Schemes for Linear Hyperbolic PDE’s 
A linear constant coefficient hyperbolic PDE has the form 



du 

dt 



a 7T - 0 

OX 



( 29 ) 



The numerical approximation of the solution at x = iAx, t = nAt is u" where Ax is the 
spatial mesh size and At is the time step. A simple first-order accurate upwinding scheme 
that takes into account the direction of wave propogation, according to the characteristic 
theory of hyperbolic PDE’s, which may be used for the solution of Eq.(29) is as follows 



Ui 






u? - Xa< 



u£i - uj 1 
u? - ujl. 



a < 0 
a > 0 



( 30 ) 



where X = At/Ax. This numerical scheme can be written in the following single equation 
form 



= uj 1 - |a[Ui n +1 - u"_ x ] + -||a|(Ui +1 - 2 Ui + uj l x ) (3ia) 



or 



ur 1 = u” - XtaMui 2 - uili) + a'luf,! - u±) ] < 31b > 

where a + = 1/2 (a + |a|), a' =1/2 (a - |a|). 
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Higher order upwind schemes may be obtained by increasing the order of accuracy of the 
stencils of Eq.(30) in the upwind direction. 

The preceding ideas may be generalized for a system of hyperbolic equations with 
constant coefficients of the following form 



du du 

dt dx 



= 0 



( 32 ) 



where U is a vector of m unknowns and A is a m x m constant coefficient matrix having 
real eigenvalues. Let X be the eigenvalue transformation diagonalizing A as X 'AX = A, 
where A is the diagonal eigenvalue matrix. Then Eq.(32) can be rearranged in diagonal 
form as 



dW 

dt 



dW 

dx 



= 0 



( 33 ) 



where W = X 'U. 

Application of the numerical schemes shown in Eq.(31) for the vector equation in 
Eq.(33) yields 



W? +1 = W? - {W? +1 - Wt i} + ||A|{ftf +1 - 2 W? - WlJ (34) 

This scheme can be written in terms of the original variables as 

C^ +1 = u? - -jA { Uj* +1 - U[ Li) + ±\A\{U&i - 2U? - Uf.J 05) 

where |A | = XAX' 1 . Similar to the one equation case the new variables A + and A' may 
be defined as A + = 1/2(A + |A|) and A' = 1/2(A - |A|). Equation (35) now takes the 
form 
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( 36 ) 



or 1 = u ; - \ {f ± . U2 - 

where 

Fj.i/2 - jMWi., * U 2 ) - lAlft;,.! - £7,)] 

?i-in = \^ u i * ~ 1^1 < u i - Wf-i)) 

The term t is referred to as the numerical flux. 

3. Conservative Schemes and Shock-Capturing Theory 
The theory of modem shock-capturing numerical schemes for non-linear hyperbolic 
conservation laws such as the Euler equations relies on the basic first-order upwinding 
methods discussed in the previous section. High order shock-capturing methods suitable 
for gas dynamic hyperbolic equations and nonconstant coefficient cases use non-linear 
differencing schemes. This significantly affects the stability of these methods because the 
stability checking can only occur locally and for the linearized versions of the non-linear 
equations. This local, linearized stability is neither necessary nor sufficient for non-linear 
problems such as shock discontinuities. Traditional methods of removing non-linear 
instability involve the introduction of linear or non-linear numerical dissipation (artificial 
viscosity or smoothing) into the difference schemes. This approach alone does not 
guarantee convergence to a physically correct solution, and higher order accurate upwind 
schemes have been developed recently to overcome the deficiencies of the artificial 
dissipation approach. 
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The limit solution of any finite-difference scheme in a conservation form consistent 
with the conservation laws satisfies the jump conditions across a discontinuity 
automatically [Ref. 14]. Weak solutions (solutions with shocks and contact 
discontinuities) of hyperbolic laws are not uniquely determined by their initial values 
except if the entropy condition, given in the previous section, is satisfied in order to pick 
out the physically relevant solution. Finite-difference approximations, however, always 
converge to a physically relevant solution when the numerical scheme used is monotone. 
If the scheme is not monotone then it must be consistent with the entropy inequality to 
assure convergence. These monotone schemes possess many desirable qualities required 
for the calculation of discontinuous solutions. 

a. Monotone Schemes 
Consider the scalar hyperbolic conservation law 

u t + au x = 0 (37) 

where a(u) = is a characteristic speed. A three-point explicit scheme in conservation 
form for the numerical integration is 

uT 1 = tij 1 - A(2£ 1/a - hf. 1/2 ) 08) 

where X = Ax/At and the numerical flux function h", 1/2 - h(u ", uj\i) is consistent with the 
conservation law in the sense that h(Uj,Uj) = f(Uj). 

The numerical scheme in Eq.(38) is of the following generalized functional form 
and it is monotone if g is a monotonically increasing function of each of its arguments. 
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uT 1 = g(ujl 1 ,u i n / ii i n +1 ) 



( 39 ) 



Monotone schemes have the attractive property of producing smooth transitions near 
discontinuities, but they are only first order accurate. Unfortunately, not all first-order 
upwind schemes are monotone. In addition, monotone, first-order upwind schemes cannot 
produce solutions for the flow field with the required accuracy using reasonable grid 
spacing. Still, higher-order shock-capturing schemes required for practical computations 
require some linear or non-linear numerical dissipation. Unfortunately, these higher order 
schemes suffer from the following deficiencies: they produce spurious oscillations when 
the solution is not smooth, they may develop instabilities at discontinuities and they may 
select nonphysical solutions. 

Two classes of modem shock-capturing schemes, the TVD and the ENO schemes, 
are more appropriate for the computation of weak solutions, especially when these 
schemes are consistent with the entropy inequality and are second-order or higher in 
smooth regions. The main distinction between ENO and TVD schemes is that the first 
can retain higher than first-order terms at extrema while the TVD schemes reduce to 
spatially first-order in these regions. However, they are more efficient in terms or 
operations count. The next section lays out the development of TVD schemes. 
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b. TVD Schemes 



Consider a family of five point difference schemes of the form 

u; f 1 + X0 (hT.ln ~ hT.l, 2 ) = u? - X(l - 6) (A£ 1/2 - h?. x/2 ) (40) 

where 

t I~1 i / J1 J1 71 73 \ 

^i + 1/2 “ ( u n-l * U 1 • u i + l • u i + 2/ 

and the scheme can be explicit or implicit depending on the value of the parameter 0. If 
0 = 0, the scheme is explicit; if 0 = 1/2, the scheme represents the trapezoid rule; and if 
0 = 1, the scheme is a backward Euler method. 

Defining 

i /2 - (i - e) h?, 1/2 * e*j ;; /2 mi 

Eq.(40) becomes 

uT 1 = uj 1 - X (F i+1/2 - h ^, 2 ) (42) 

which is of the form 

= L 2 u D 

with L, and Lj defined as the following difference operators 

L a u = u 1 + A.6(i2 i+;L/2 - h d _ 1/2 ) 

L 2 u = u ± + X( 1-0) (i2 i+1/2 - Ai_ 1/2 ) 

The total variation of the mesh u" is defined as 
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( 43 ) 



TV(u n ) = E|ujU - ui 1 ! = S|A i+1/2 u 
The numerical scheme is TVD if the following inequality holds 



7V(u n+1 ) <; TV( u n ) 

There are four design principles used in the construction of high-resolution TVD 
schemes: 

• The flux-corrected transport of Boris and Book [Ref. 15]. This is a two step hybrid 
scheme consisting of a combined first and second-order scheme. It computes an 
update form of a first-order scheme and then filters the second-order corrections by 
the use of flux emitters to prevent occurence of new extrema. 

• A second-order extension of Gudunov’s scheme by van Leer [Ref. 16] and Colella 
and Woodward [Ref. 17]. This scheme is based on the observation that one can 
obtain second-order accuracy in Gudunov’s scheme by replacing the piecewise- 
constant initial data of the Reimann problem with piecewise-linear or parabolic 
initial data. The slope of the initial data is chosen so that no spurious oscillations 
can occur. 

• The modified-flux TVD scheme of Harden [Ref. 18]. This is a technique to design 
a second-order accurate TVD scheme by starting with a first-order TVD scheme and 
applying it to a modified flux. The modified flux is chosen so that the scheme is 
second-order in smooth regions and first-order at extrema. 

• The numerical fluctuation approach of Roe [Ref. 19] depending on an average 
function. The average function is constructed in such a way that spurious 
oscillations will not occur by the use of flux limiters. 

Most of the above methods can be viewed as three-point central-difference schemes 
with an appropriate numerical dissipation term which automatically controls the amount 
of numerical dissipation. This control of the dissipation is automatic unlike the control 
of the numerical dissipation used by central-differencing techniques used in linear theory. 
The basic idea of the above design principles is to construct a higher-order scheme which 
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has properties similar to a first-order TVD scheme so that spurious oscillations cannot be 
generated. In order to achieve that limiting procedure, a flux limiter is used, which 
imposes constraints on the gradients of the dependent variables or the flux function. 

The derivation of the higher-order schemes with TVD properties can be achieved 
by either following the so-called MUSCL (Monotonic Upstream Scheme for Conservation 
Laws) approach of van Leer [Ref. 16], Harden [Ref. 18], Roe [Ref. 19] and Osher and 
Chakravarthy [Ref. 20]; and the non-MUSCL approaches of Roe [Ref. 21] and Yee 
[Ref. 22]. 

4. Development of Flow Simulation Code 

The numerical scheme used for the simulation was developed in Ref. 23. A 
summary of the numerical scheme is given next. 

The general form of the conservation law form for an inertial reference frame is 



The main motivation of the conservation law form is to capture flow discontinuities. 
Finite volume discetization is obtained by replacing the surface integrals by sums over 
the faces of a cell. The differential form is obtained by applying Eq.(44) to a cell in 
physical space. In order to facilitate the numerical computation for complex geometries, 
the physical space cell (x,z) is the image of the computational domain cell (£,£) resulting 
from the generalized coordinate transformation 




( 44 ) 



v 



s 



V 



X = x(££,t) , z = z(££,t) 



(45) 
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The transformed coordinate system (£,Q is an evenly spaced rectangular domain 
(orthogonal) having area of unity. This allows all numerical differentiations to be 
performed with unweighted differencing formulas. In addition, the conservation law form 
represents the thin-layer approximation. Terms of order 1/Re 1/2 and smaller are neglected, 
and the viscous terms containing derivatives in the streamwise direction are dropped. 
These viscous terms are smaller in magnitude than the viscous terms with derivatives 
normal to the body surface. The thin-layer approximation allows a reduction in 
computational time, since the streamwise viscous gradients are not usually resolved 
adequately anyway. Also, the thin-layer approximation allows computation of separated 
flows and flows with normal pressure gradients. The conventional boundary layer 
approximation does not allow these computations. 

The differential form of the two-dimensional conservation law in generalized 
coordinates is as follows: 



d t Q + d + d„F = Re-'d^S ( 46 ) 

where the transformed vectors are 





'P' 




' pu ' 




f pW > 


~ 1 


pu 


A 1 


puU + 


~ I 


P UW + CxP 


Q = — 
J 


pw 


• ? 


pvU + IjD 


' F -7r 


pww + C 2 P 








K (e + p)u - l t p ) 




^(e + p) - C cPj 



95 



/ 



0 > 

^jU c + ( M- / 3 ) ni 2 C x 
+ ( M- / 3 ) /n 2 

p/n^ + (p/3 )/n 2 + (C x u + C y v + C z w) ; 

1 lerc 

/n a = Cx + C* 

= C x u c + Oc < 48 > 

/n, = (u 2 + w 2 ) / 2 + KPr' 1 ( °_ F - ) 

dC 

and U and W arc the contravariant velocity components given by 

U = Ul x + w£ z + ( 49 ) 

W = uC x + + c t 

In the above equations, all geometric dimensions arc normalized with the chord length; 
p is the density normalized with the freestream density p_; u and w arc the Cartesian 
velocity components normalized with the speed of sound; e is the total energy per unit 
volume normalized with p_a‘_; Pr is the Prandtl number; and K is the thermal 
conductivity. The pressure is related to the density and total energy through the equation 
of state for an ideal gas, p = (y - l)[c - p(u 2 + w 2 )/2]. 

Spatial schemes for the numerical solution of the inviscid part of Eq. (46) are of the 
following semidisercte form 
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dQ _ _ E i , Jc+1/2 i , k-1/2 _ Jc+1/2 ^i.k- 1/2 ,^ni 

dt AC AC 

where H and G are the fluxes defined as 

E 1*i/ 2, k ~ \ t ( E i.k ~ E l+l,k) ~ A |A‘li + i/2,JcJ 

with 

A \ E \i*i/2.k ~ \A\(Ql+i.k ~ 0 ltk ) 

1^ I j+ 1 / 2 , k = MO) 'Q = 0(0 Ul . k - Qi, k ) 

Flux difference splitting can be accomplished by enforcing interaction of neighboring 
states using Roe’s approximate Reimann solver, where the interaction is enforced through 
neighboring cell interfaces. The local time linearization of the difference vectors AE* for 
left and right states Ql and Qr, respectively, using Roe’s method is obtained as follows 

A E* = A{Q) t LQ , A* = RA^R -1 < 51 > 

where AQ = Qr - Ql, A = (3 £/3£) are the flux Jacobian matrices. A 1 is the eigenvalue 
diagonal matrix and RJR ' 1 are the left and right eigenvector matrices. The intermediate 
state 0 is calculated using the Roe variables as 
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p = /pIp^ 

- = “l/P ~l + “sy/P^ 

\/Pl + y/f*R 

w = - Wr ^~ r ( 52 ) 

yj Pi + 

_ ^l/Pl + ^r^Pr 
\[P~L + \fP~R 

The numerical integration scheme is performed using a fully-upwind, factored, implicit, 
third order accurate numerical scheme. The upwinding is perfomed with flux-difference 
splitting using Roe’s scheme. For complex flow fields first order accurate schemes are 
too diffusive and they do not produce solutions with reasonable grid spacings. Here a 
higher order scheme is used. The numerical scheme is given below: 

[I- + V1&.JC + 4s;. J =>] p x 1531 

[I * + Af. c;.* - se-'SfMi,*)]" x (Oft - Sf. k ) = 

- [@l.k ~ &l,k + ^i.1/2,* “ + fy(^i,Jc,l/2 “ ^l.k-l/2^ p 

-Re “ ^l,k-l/2^ P ^ 

= Ax/A^, and A, V and 5 are the forward, backward and the central difference 
operators, respectively. The quantities S itU2 , k , P i>A , l/2 , are numerical fluxes. 
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The accuracy of the numerical solution is improved by subiterations to convergence 



within each time step iteration. The approximation to O'" 1 at each subiteration is Q p . 
When p > 2, during a given level of subiteration to convergence, 0 P = 0"* 1 ; but when p 
= 1 and no subiterations are performed, then 0 P = 0\ and 0 P ' 1 = O'" 1 . By subiterating to 
convergence, linearization and factorization errors can be eliminated during the iteration 
process, because the left hand side of Eq.(53) can be driven to zero. The inviscid fluxes £ 
and F are evaluated using the Roe upwinding. The numerical fluxes for a third order 
accurate scheme are given by 

f i.k = h(Q itlik .Qj, k ) ~ [dfi. 3 / 2 . Jt + 2df; +1/2fJt ] (54) 

+ [2cffj +1 / 2 ,^ + df±_ 1 /2 l k ] 

or in terms in expanded form 

1/2, k ~ ~2 [EiQi'kt ti+i/2. k ) ~ E(Qi-\. k > ^J+ 1 / 2 ,*)^ (55) 



c (Oi tk i Pi+i,ic' + 2AE (Qi+ 2 , k > Pi+i,lc' ^i+3/2, 



+ [2A E (Qi +lik , Qi ik i £i-i/2,/c) + (Ql, k ' Qi+I, k > £i-i/2,Je) ] 



Using the notation of Eq.(54), h corresponds to the first order accurate flux, and the df 
and df terms represent the correction terms added to the first order accurate fluxes. 
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Higher order accurate shock-capturing schemes have some limitations: they may select 
a nonphysical solution, they produce spurious oscillations and they may develop nonlinear 
instability in nonsmooth and discontinuous flow regions. More appropriate shock- 
capturing schemes suitable for the computation of weak solutions are the TVD schemes, 
described earlier. For this simulation, the Osher-Chakravarthy TVD scheme [Ref. 20] was 
implemented. This scheme uses flux limiters to impose constaints on the gradient of the 
fluxes. The flux-limited values are computed as follows: 

3 f 1*3/2, k = minmod[dfi+2/ 2l ki 1*1/2 . *] 

3 fi +1/2, k minmod [cf.fj+i/2, $df±+2/2,k\ ($6) 

1*1/2, k ~ minmod [dfj +1 / 2 ,k 1 1-1/2, 

1-1/2, k = minmod[dfi- 1 /2,kr$dfi+ 1 /2,k] 

where the minmod operator is defined by 

minmod[x,y ] = sign(x) x max[0,min{ /x /,ysign(x)}J (57) 

and the comparison parameter 0 is in the range 1 s p s 4. Substitution of the flux 
limited values given by Eq.(56) into Eq.(54) makes the entire scheme TVD. 

The flux differences AE are evaluated by the Roe scheme described earlier. The 
viscous fluxes s i k , Ui are computed using central differences as shown below 
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( 58 ) 



^l,k*l/2 [Qi, Jc+1/2 ' (@0 l,k*l/2' ^l.k+1/2] 

@1, k+l/2 ~ (Ql.k ~ @1+1, k^ 

{QO 

1, k+l/2 “ @l,k+l @l,k 

The eigenvalues of the flux Jacobian matrices A ,B are evaluated as suggested in Ref. 
24 to prevent expansion shocks. The terms A ,C of the left hand side operators of 
Eq.(53) have been linearized using the Steger-Warming flux vector splitting. A validation 
of the viscous flow algorithm for a transonic airfoil is shown in Fig. D2. 
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d3- 



Figure D1 . Comparison of Turbulence Model Predictions [Ref. 23] 
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APPENDIX E 



GRAPE 



The following is a summary of the grid-generating algorithm GRAPE. Reference 
25 provides a detailed explanation of the algorithm theory and Ref. 26 is a detailed user’s 
manual. GRAPE is an algorithm that transforms x and y coordinates in real space to £ 
and r| coordinates in computational space. GRAPE stands for GRids about Airfoils using 
Poisson’s Equation, rj = 0 forms the inner boundary of the grid, and r| = forms the 
outer boundary. For the C-grids, £ = 0 is located at the rear boundary, and £ moves 
forward clockwise around the airfoil until it reaches at the rear boundary again. The 
mapping functions are required to satisfy the following: 



+ % y = p ( la ) 

+*lyy=Q (lt>) 



The following relations are useful in transforming the equations from real space to 
computational space: 



where 



II 

V- 


(2a) 


£y = -V J 


(2b) 


n* = -y&/J 


(2c) 


riy = x ? /J 


(2d) 


j = 


(2c) 
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Applying Eqns. (2) to (1) yields the following transformed Poisson Equations: 



<JXf£ - 2|ix^ + = -J 2 (Px^ + Qx„) 

- 2|iy^ + ty^ = -J 2 (Py^ +Qy„) 



(3b) 



(3a) 



where 



o = x„ 2 + y„ 2 



(3c) 



P = + y^ 



(3d) 



t = x ^ 2 + y ^ 2 



(3e) 



Solving Eqn. (3) for a particular choice of inhomogeneous terms P and Q (forcing 
functions) and a particular set of boundary conditions causes the grid to be generated. 

There is latitude that exists with the choices of P and Q. Different choices of P and 
Q result in different grids. In an effort to choose P and Q to allow reasonable effort in 
the computation, P and Q are redefined in terms of four new variables. In addition to 
geometric constraints and the original Poisson equations, six equations with six unknowns 
are left. These equations can be solved in a straightforward iterative scheme. 

The following geometric constraints are imposed on the grid: 

1. Spacing along £ = constant lines between the body at tj = 0 and next grid node 
T| = 1 is specified by the user 

2. The angle of intersection between the body and £ = constant lines is specified by 
the user 

3. Same geometric constraint as 1) but applied to the outer boundary 

4. Same geometric constraint as 2) but applied to the outer boundary 

These geometric constraints are embedded in the P and Q terms. 
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The program uses a method called Coarse-Fine Sequencing (CFS). This method 
greatly accelerates the numerical convergence of the program. The coarse solution solves 
the grid at every third point in the a direction and every third point in the 6 direction. 
The fine solution uses a cubic spline fit to find the rest of the points. 

Reference 27 details changes made to the original GRAPE code to allow the 
generation of periodic cascade C-grids for turbomachinery flow. In addition, the outer 
and inner boundary routines were modified. All the features of the original GRAPE code 
were retained. 
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APPENDIX F 



GRAPE INPUT FILES 



$GRID1 

JMAX=199 , KMAX=64 , NTETYP=3 , NAIRF=5 , JAIRF=3 16 , MIBDST=3 , 

NDS=2 , DSI = 0 .001 , JTEBOT=l 3 , JTETOP=18 7 , XLE=0 , XTE = 1 , NOBSHP=4 ,NOBDST=l, 
XLEFT=-0. 3 , 

XRIGHT=1 . 5 , YBOTOM=-0 .276, YTOP=0 .276, RCORN=0 .20, ALAMF=1 . 0 , ALAMR=0 . 0, 
N0RDA=4 , 1 , 

MAXITA=10,10, JPRT=*-1,NOUT=0 

$END 

$GRID2 

XTFRAC=1 . 0 , ROTA1JG=0 . 0 , WAKEP=1 . 0 

$E11D 

$GRID3 

AIRFX « 6. 00000, 5. 99948 , 5.99891, 5. 99815, 5. 99667 , 
5.99481,5.99264,5.99106,5.98939,5.98766, 

5.98500, 

5.98200,5.97870,5.97507,5.97108,5.96668, 
5.96185,5.95654,5.95069,5.94426,5.93719 , 
5.92941,5.92085,5.91143,5.90108 ,5.88968, 
5.87715,5.86337,5.84820,5.83152,5.81318, 
5.79299,5.77079,5.74637,5.71951,5.68996, 
5.65745,5.62170,5.58237,5.53911,5.49152, 

5.43917 ,5.38159,5.31825,5.25556,5.19286, 
5.13016,5.06747,5.00477, 4.94208,4.87938, 
4.81668,4.75399,4.69129,4.62860,4.56590, 
4.50320,4.44051,4.37781,4.31512,4.25242, 
4.18973,4.12703,4.06433,4.00164,3.93894, 
3.87625,3.81355,3.75085,3.68816,3.62546, 

3 . 56277,3.50007,3 .43738,3.37468,3.31198, 
3.24929,3.18659,3.12390,3.06120,2.99850, 
2.93581,2.87311,2.81042,2.74772,2. 68503, 
2.62233,2.55963,2.49694,2.43424,2.37155, 
2.30885,2.24616,2.18346,2.12076,2.05807, 
1.99537,1.93268,1.86998,1.80728,1.74459, 
1.68189,1.61920,1.55650,1.49381,1.43111, 

1.36841,1.30572, 1.24302,1.18033,1.11763, 
1.05493,0.99224,0.92954 ,0.86685,0.80415, 

0.74146, 

0.67876,0.61541,0.55783,0.50548,0.45789, 

0.41460,0.37530,0.33955,0.30704,0.27749, 

0.25063,0.22621,0.20401,0.18383,0.16548, 

0.14880,0.13363,0.11985,0.10732,0.09593, 

0.08557,0.07615,0.06759,0.05981,0.05274, 

0.04631,0.04046,0.03515,0.03032,0.02592, 

0.02193,0.01830,0.01500, 

0.01369,0.01112,0.00987,0.00750,0.00536, 

0.00351,0.00201,0.00091,0.00023, 

0.00000,0.00023,0.00091,0. 00201, 0.00351, 
0.00536,0.00750,0.00987,0.01112,0.01369, 
0.01699,0.02062,0.02461,0.02901,0.03384, 

0.03915,0.04500,0 . 05143,0.05850,0.06628, 
0.07483,0.08426,0.09462 ,0.10601,0.11854 , 
0.13232,0.14749,0.16417,0.18252,0.20270, 
0.22490,0.24932,0.27887,0.31138,0.34713, 
0.38646,0.42972,0.47731,0.52966,0.58724 , 

0.65058 ,0.71308,0.77608,0.83908,0.90208, 

0.96508, 1.02508 ,1.09108,1.15108,1.21408, 
1.27708,1.34008,1.40308,1.46608,1. 52908, 

1.59208 ,1.65208,1.71508,1.77808, 1.84108, 
1.90408,1.96708,2.03008,2.09308,2.15608, 



Figure FI. GRAPE Input for Channel C-grid 
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2.21908,2.28208,2.34208,2.40508,2.46808, 

2.53108,2.59408,2.65708,2.72008,2.78308, 

2.85000,2.92118,2.99378,3.05371,3.11665, 

3.17956,3.24256,3.30553,3.36851,3.43149, 

3.49449,3.55748,3.61749,3.68049,3.74349, 

3.80648,3.86948,3.93247,3.99546,4.05844, 

4.12141,4.18137,4 .24431,4.30725,4.37108, 

4.43309,4.49598,4.55886,4.62171,4.68455, 

4.74736,4.81015,4.87292,4.93566,4.99837, 

5.06105,5.12371,5.18633, 

5.24892,5.31849,5.38183, 

5.43945,5.49176,5.53935,5.58262,5.62195, 

5.65770,5.68899,5.71975,5.74662,5.77104, 

5.79324,5.81350,5.83177,5.84845,5.86361, 

5.87740,5.88993,5.90132,5.91168,5.92109, 

5.92965,5.93743,5.94451,5.95094,5.95678, 

5.96210,5.96741,5.97224,5.97664,5.98063, 

5.98426, 

5.98756,5.98930,5.99097,5.99256,5.99474, 
5.99603,5.99810,5.99888,5 .99945,6.00000, 
AIRFY = 0.01500,0.01109,0.00940,0.00779,0.00557, 
0.00365,0.00209,0.00128,0.00066,0.00024 , 
0.00000, 

0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0.00000,0.00000,0.00000,0.00000,0 . 00000, 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0.00000,0.00000,0.00000,0.00000,0.00000, 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0.00000,0.00000,0.00000,0.00000,0.00000, 
0.00000,0.00000,0.00000,0.00000,0.00000, 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0.00000,0.00000,0.00000,0.00000,0.00000, 
0.00000,0.00000,0.00000,0.00000,0.00000, 
0.00000,0.00000,0.00000,0.00000,0.00000, 
0.00000,0.00000,0.00000,0.00000,0.00000, 
0.00000,0.00000,0.00000,0.00000,0.00000, 
0.00000, 

0.00000,0.00000,0.00000,0.00000,0.00000, 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0.00000,0.00000,0.00000,0.00000,0.00000, 
0.00000,0.00000,0.00000,0.00000,0.00000, 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0.00000,0.00000,0.00000, 

0.00006,0.00051,0.00091,0.00201,0.00351, 
0.00536,0.00750,0.00987,0.01240,0.01500, 
0.01761,0.02013,0.02250,0.02464,0.02640, 
0.02910,0.02940,0.02994 , 

0.03015,0.03034,0.03061,0.03089,0.03118, 

0.03150,0.03186,0.03225,0.03272,0.03316, 

0.03364,0.03426,0.03493,0.03559,0.03639, 

0.03720,0.03813,0.03915,0.04025,0.04150, 



Figure FI (cont) . GRAPE Input for Channel C-grid 
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0 . 04282 , 0 . 04436 , 0 . 04612 , 0 . 04814 , 0 . 05034 , 

0 . 05272 , 0 . 05539 , 0 . 05830 , 0 . 06153 , 0 . 06502 , 

0 . 06887 , 0 . 07273 , 0 . 07658 , 0 . 08043 , 0 . 08428 , 

0 . 08814 , 0 . 09181 , 0 . 09584 , 0 . 09951 , 0 . 10337 , 

0 . 10722 , 0 . 11107 , 0 . 11493 , 0 . 11878 , 0 . 12263 , 

0 . 12649 , 0 . 13016 , 0 . 13401 , 0 . 13786 , 0 . 14172 , 

0 . 14557 , 0 . 14942 , 0 . 15328 , 0 . 15713 , 0 . 16098 , 

0 . 16484 , 0 . 16869 , 0 . 17236 , 0 . 17621 , 0 . 18007 , 

0 . 18392 , 0 . 18777 , 0 . 19163 , 0 . 19548 , 0 . 19933 , 

0 . 20342 , 0 . 20763 , 0 . 21145 , 0 . 21434 , 0 . 21709 , 

0 . 21955 , 0 . 22172 , 0 . 22359 , 0 . 22517 , 0 . 22645 , 

0 . 22745 , 0 . 22814 , 0 . 22854 , 0 . 22866 , 0 . 22850 , 

0 . 22804 , 0 . 22728 , 0 . 22623 , 0 . 22489 , 0 . 22326 , 

0 . 22133 , 0 . 21923 , 0 . 21673 , 0 . 21394 , 0 . 21085 , 

0 . 20747 , 0 . 20380 , 0 . 19984 , 0 . 19558 , 0 . 19103 , 

0 . 18619 , 0 . 18106 , 0 . 17563 , 0 . 16991 , 0 . 16390 , 

0 . 15760 , 0 . 15101 , 0 . 14412 , 

0 . 13695 , 0 . 12859 , 0 . 12073 , 

0 . 11328 , 0 . 10630 , 0 . 09977 , 0 . 09367 , 0 . 08802 , 

0 . 08278 , 0 . 07811 , 0 . 07345 , 0 . 06931 , 0 . 06551 , 

0 . 06201 , 0 . 05879 , 0 . 05585 , 0 . 05316 , 0 . 05064 , 

0 . 04841 , 0 . 04634 , 0 . 04444 , 0 . 04271 , 0 . 04113 , 

0 . 03969 , 0 . 03781 , 0 . 03717 , 0 . 03608 , 0 . 03508 , 

0 . 03417 , 0 . 03326 , 0 . 03242 , 0 . 03167 , 0 . 03098 , 

0 . 03035 , 

0 . 02978 , 0 . 02937 , 0 . 02876 , 0 . 02796 , 0 . 02641 , 

0 . 02451 , 0 . 02230 , 0 . 02069 , 0 . 01901 , 0 . 01500 , 

$END 



Figure FI (cont) . GRAPE Input for Channel C-grid 
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$GRID1 




JMAX=169 , KMAX=4 0 , NTETYP= 


3, HAIRF*=5 , HIBDST=7, MOBSHP=7, 


JAIRF=316 , JTEBOT= 25, JTETOP=145, NORDA=0,3, MAXITA= 0,50, NOUT=4 , 


DSI=. 0025 , XLE=0 . 0 , XTE=0. 


61786, 


XLEFT=-0 .3089, XRIGHT=1 . 2 35 , RCORR* 0.0333, 


$END 




$GRID2 




N0BCAS=0, NLE=2 0 , NTE= 10, 


DSRA= 0.5, 


DSLE= 0.0005, DSTE= 0.0010 


i, PITCH= 0.50, 


YSCL=1 . 0 , XTFRAC=0 . 8 , ROTANG= -51.84, 


WAKEP= 0.8, 




$EIJD 




$GRID3 




AIRFX= 




6.00000, 5.99948, 


5.99891, 5.99815, 5.99667 


5.99481, 5.99264, 


5.99106, 5.98939, 5.98766, 


5.98500, 




5.98200, 5.97870, 


5.97507,5.97108,5.96668, 


5.96185, 5.95654, 


5.95069,5.94426,5.93719, 


5.92941, 5.92085, 


5.91143,5.90108,5.88968, 


5.87715, 5.86337, 


5.84820,5.83152,5.81318, 


5.79299, 5.77079, 


5.74637,5.71951,5.68996, 


5.65745, 5.62170, 


5.58237,5.53911,5.49152, 


5.43917, 5.38159, 


5.31825,5.25556,5.19286, 


5.13016, 5.06747, 


5.00477,4.94208,4.87938, 


4.81668, 4.75399, 


4.69129,4.62860,4.56590, 


4.50320, 4.44051, 


4.37781,4.31512,4.25242, 


4.18973, 4.12703, 


4.06433,4.00164,3.93894, 


3.87625, 3.81355, 


3.75085,3.68816,3.62546, 


3.56277, 3.50007, 


3.43738,3.37468,3.31198, 


3.24929, 3.18659, 


3 . 12 390,3.0612 0,2.99850, 


2.93581, 2.87311, 


2.81042,2.74772,2.68503, 


2.62233, 2.55963, 


2.49694,2.43424,2.37155, 


2.30885, 2.24616, 


2 . 18 34 6,2.12 076,2.058 07, 


1.99537, 1.93268, 


1.86998,1.80728,1.74459, 


1.68189, 1.61920, 


1.55650,1.49381,1.43111, 


1.36841, 1.30572, 


1.24302,1.18033,1.11763, 


1.05493, 0.99224, 


0.92954,0.86685,0.80415, 


0.74146, 




0.67876, 0.61541, 


0.55783,0.50548,0.45789, 


0.41460, 0.37530, 


0.33955,0.30704,0.27749, 


0.25063, 0.22621, 


0.20401,0.18383,0.16548, 


0.14880, 0.13363, 


0.11985,0.10732,0.09593, 


0.08557, 0.07615, 


0.06759,0.059 81,0.0527 4 , 


0.04631, 0.04046, 


0.03515,0.03032,0.02592, 


0.02193, 0.01830, 


0.01500, 


0.01369, 0.01112, 


0.00987,0.00750,0.00536, 


0.00351, 0.00201, 


0.00091,0.00023, 


0.00000, 0.00023, 


0.00091,0.00201,0.00351, 


0.00536, 0.00750, 


0.00987,0.01112,0.01369, 


0.01699, 0.02062, 


0.02461,0.02901,0.03384, 


0.03915, 0.04500, 


0.05143,0.05850,0.06628, 


0.07483, 0.08426, 


0.09 4 62,0.10601,0.11854 , 


0.13232, 0.14749, 


0.16417,0.18252,0.20270, 


0.22490, 0.24932, 


0.27887,0.31138,0.34713, 


0.38646, 0.42972, 


0.47731,0.52966,0.58724, 


0.65058, 0.71308, 


0.77608,0.83908,0.90208, 


0.96508, 1.02508, 


1.09108,1.15108,1.21408, 


1.27708, 1.34008, 


1.40308,1.46608,1.52908, 


1.59208, 1.65208, 


1.71508,1.77808,1.84108, 



Figure F2. GRAPE Input for Cascade C-grid 
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AIRFY= 



1 . 90408 , 
2 .21908 , 
2 . 53108 , 
2 . 85000 , 
3 . 17956 , 
3 . 49449 , 
3 . 80648 , 
4 . 12141 , 
4 . 43309 , 
4 . 74736 , 
5 . 06105 , 
5 . 24892 , 
5 . 43945 , 
5 . 65770 , 
5.79324 , 
5 . 87740 , 
5 . 92965 , 
5 . 96210 , 
5 . 98426 , 
5 . 98756 , 
5 . 99603 , 



1 . 96708 , 
2 . 28208 , 
2 . 59408 , 
2 . 92118 , 
3 . 24256 , 
3 . 55748 , 
3 . 86948 , 
4 . 18137 , 
4 . 49598 , 
4 . 81015 , 
5 . 12371 , 
5 . 31849 , 
5 . 49176 , 
5 . 68899 , 
5 . 81350 , 
5 . 88993 , 
5 . 93743 , 
5 . 96741 , 

5 . 98930 , 

5 . 99810 , 



2 . 03008 
2 .34208 
2 . 65708 
2 . 99378 
3 . 30553 
3 .61749 
3 . 93247 
4 .24431 
4 . 55886 
4 .87292 
5.18633 
5.38183 
5.53935 
5.71975 
5.83177 
5.90132 
5.94451 
5.97224 

5.99097 

5.99888 



, 2 . 09308 
, 2 .40508 
, 2 .72008 
, 3 . 05371 
, 3 . 36851 
, 3 . 68049 
, 3 .99546 
, 4.30725 
, 4 . 62171 
, 4.93566 



, 5.58262 
, 5.74662 
, 5. 84845 
, 5.91168 
, 5.95094 
, 5.97664 

, 5.99256 

, 5.99945 



, 2 . 15608 , 
, 2 . 46808 , 
, 2 . 78308 , 
, 3 . 11665 , 
, 3.43149 , 
, 3 . 74349 , 
, 4 . 05844 , 
, 4 . 37108 , 
, 4 . 68455 , 
, 4 . 99837 , 



, 5 . 62195 , 
, 5 . 77104 , 
, 5 . 86361 , 
, 5 . 92109 , 
, 5 . 95678 , 
, 5 . 98063 , 

, 5 . 99474 , 

, 6 . 00000 , 



0 . 01500 , 0 . 01109 , 0 . 00940 , 0 . 00779 , 0 . 00557 , 

0 . 00365 , 0 . 00209 , 0 . 00128 , 0 . 00066 , 0 . 00024 , 

0 . 00000 , 

0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 

0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 
0 . 00000 , 

0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 

0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 

0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 

0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 

0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 

0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 0 . 00000 , 

0 . 00000 , 0 . 00000 , 0 . 00000 , 

0 . 00006 , 0 . 00051 , 0 . 00091 , 0 . 00201 , 0 . 00351 , 
0 . 00536 , 0 . 00750 , 0 . 00987 , 0 . 01240 , 0 . 01500 , 
0 . 01761 , 0 . 02013 , 0 . 02250 , 0 . 02464 , 0 . 02640,0 
0 . 02910 , 0 . 02940 , 0 . 02994 , 

0 . 03015 , 0 . 03034 , 0 . 03061 , 0 . 03089 , 0 . 03118 , 

0 . 03150 , 0 . 03186 , 0 . 03225 , 0 . 03272 , 0 . 03316 , 
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. 02799 , 
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0 . 03364 , 0 . 03426 , 0 . 03493 , 0 . 03559 , 0 . 03639 , 

0 . 03720 , 0 . 03813 , 0 . 03915 , 0 . 04025 , 0 . 04150 , 

0 . 04282 , 0 . 04436 , 0 . 04612 , 0 . 04814 , 0 . 05034 , 

0 . 05272 , 0 . 05539 , 0 . 05830 , 0 . 06153 , 0 . 06502 , 

0 . 06887 , 0 . 07273 , 0 . 07658 , 0 . 08043 , 0 . 08428 , 

0 . 08814 , 0 . 09181 , 0 . 09584 , 0 . 09951 , 0 . 10337 , 

0 . 10722 , 0 . 11107 , 0 . 11493 , 0 . 11878 , 0 . 12263 , 

0 . 12649 , 0 . 13016 , 0 . 13401 , 0 . 13786 , 0 . 14172 , 

0 . 14557 , 0 . 14942 , 0 . 15328 , 0 . 15713 , 0 . 16098 , 

0 . 16484 , 0 . 16869 , 0 . 17236 , 0 . 17621 , 0 . 18007 , 

0 . 18392 , 0 . 18777 , 0 . 19163 , 0 . 19548 , 0 . 19933 , 

0 . 20342 , 0 . 20763 , 0 . 21145 , 0 . 21434 , 0 . 21709 , 

0 . 21955 , 0 . 22172 , 0 . 22359 , 0 . 22517 , 0 . 22645 , 

0 . 22745 , 0 . 22814 , 0 . 22854 , 0 . 22866 , 0 . 22850 , 

0 . 22804 , 0 . 22728 , 0 . 22623 , 0 . 22489 , 0 . 22326 , 

0 . 22133 , 0 . 21923 , 0 . 21673 , 0 . 21394 , 0 . 21085 , 

0 . 20747 , 0 . 20380 , 0 . 19984 , 0 . 19558 , 0 . 19103 , 

0 . 18619 , 0 . 18106 , 0 . 17563 , 0 . 16991 , 0 . 16390 , 

0 . 15760 , 0 . 15101 , 0 . 14412 , 

0 . 13695 , 0 . 12859 , 0 . 12073 , 

0 . 11328 , 0 . 10630 , 0 . 09977 , 0 . 09367 , 0 . 08802 , 

0 . 08278 , 0 . 07811 , 0 . 07345 , 0 . 06931 , 0 . 06551 , 

0 . 06201 , 0 . 05879 , 0 . 05585 , 0 . 05316 , 0 . 05064 , 

0 . 04841 , 0 . 04634 , 0 . 04444 , 0 . 04271 , 0 . 04113 , 

0 . 03969 , 0 . 03781 , 0 . 03717 , 0 . 03608 , 0 . 03508 , 

0 . 03417 , 0 . 03326 , 0 . 03242 , 0 . 03167 , 0 . 03098 , 

0 . 03035 , 

0 . 02978 , 0 . 02937 , 0 . 02876 , 0 . 02796 , 0 . 02641 , 

0 . 02451 , 0 . 02230 , 0 . 02069 , 0 . 01901 , 0 . 01500 , 

SEND 
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